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ABSTRACT 



Context. Only a small number of high mass stars (> 30 M Q ) have fundamental parameters (i.e. masses and radii) measured with high 
enough accuracy from eclipsing binaries to constrain formation and evolutionary models of massive stars. 

Aims. This work aims to increase this limited sample, by studying the four massive eclipsing binary candidates discovered by Bonanos 
in the young massive cluster Westerlund 1. 

Methods. We present new follow-up echelle spectroscopy of these binaries and models of their light and radial velocity curves. 
Results. We obtain fundamental parameters for the eight component stars, finding masses that span a range of 10 - 40 M Q , and con- 
tributing accurate fundamental parameters for one additional very massive star, the 33 M Q component of W13. WR77o is found to 
have a ~ 40 M Q companion, which provides a second dynamical constraint on the mass of the progenitor of the magnetar known in 
the cluster. We also use W13 to estimate the first, direct, eclipsing binaiy distance to Westerlund 1 and therefore the magnetar and 
find it to be at 3.7 ± 0.6 kpc. 

Conclusions. Our results confirm previous evidence for a high mass for the progenitor of the magnetar. In addition, the availabil- 
ity of eclipsing binaries with accurate parameters opens the way for direct, independent, high precision eclipsing binary distance 
measurements to Westerlund 1 . 

Key words, stars: fundamental parameters, stars: massive, binaries: eclipsing, stars: early-type, stars: Wolf-Rayet, open clusters and 
associations: individual (Westerlund 1) 



1. Introduction 

Improving stellar structure, formation and evolution models of 
massive stars is limited by the accuracy with which their fun- 
dame ntal parameters (masses and radii) are known. lAnderseril 
d 1 99 lb reviewed the im portance of accura te parameters from 
eclipsing binaries, while iTorres et"aD (1201 Ol) revisited the sub- 
ject, reviewing the advances made possible by the continuous 
improvements in observational accuracy. They also compiled a 
catalog of detached systems with measurements accurate to +3% 
from the literature, which consists of 95 systems and extends the 
mass range from O.2M to 27 M . The earliest -type star in the 
sample was an 07 star with a mass of 27.3 Mn. iBonanosi (f2009) 
found that the less stringent accuracy criterion of 10% is only 
fulfilled by 15 eclipsing binaries with masses above 30 M , this 
number having only increas ed to 17 (with th e components of 
Cyg OB2-17, measured by IStroud et al .1 120101) in the last three 
years. This astounding lack of accurate measurements of the fun- 
damental parameters of the highest mass stars, combined with 
the large parameter space of metallicity and evolutionary stage 
that they span, which is further complicated by multiplicity and 
in particular binary evolution, serves as the motivation for this 
work. 

To improve the situation, a systematic survey of the most 
massive stars in the Local Group has been undertaken (see 
iBonanoj 120101 for an overview). As part of this broader sur- 
vev. lBonanosi (2007) conducted the first variability survey of the 



young massive "super star cluster" Westerlund 1 (IClark et al.1 
2005) in search of massive eclipsing binary candidates, as it 
is one of the most massive youn g clusters known in the Local 
Group, with an age of 4.5-5 Myr (ICrowther e t al. 2006). Wd 1 is 
of great interest, as it contains an assortment of rare evolved mas- 
sive stars, such as blue, ye llow and red supergiants, including 
a rare supergiant B[e] star (IClark et al. 2005t iNeeueruela et al.l 
1201 0l) . 24 confirmed W olf-Rayet stars (ICrowther et al.ll2006l) . a 
lumin ous blue variable (IClark & Negueruelall2004l) . and a mag- 
netar (Mu no et"a"fll2006l) . thereby making possible the determi- 
nation of fundamental parameters of massive stars at different 
evolutionary stages. The photometric variability search resulted 
in four massiv e eclipsing bin ary candidates (Wdeb, W13, W36, 
and WR77o: lBonanosl2007l) . which are the subject of this work. 

In the following sections, we present follow-up spectroscopy 
and modeling of these four binaries, and the resulting component 
parameters. We present the observations in Section 2, the radial 
velocity measurements in Sect. 3, the models for each eclipsing 
binary in Sect. 4, the comparison with evolutionary models in 
Sect. 5, an estimation of the distance in Sect. 6, and a discussion 
of our results and conclusions in Sect. 7. 
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2. Observations and data reduction 

Multi-epoch spectroscopy of the four eclipsing binaries was ob- 
tained during seven nights in 2 007-2008 and one nig ht in 2011 
using the MIKE spectrograph dBernstein et alj|2~003l) at the 6.5 
meter Magellan (Clay) telescope at Las Campanas Observatory, 
Chile. This camera uses a 2048 x 4096 SITe CCD with a 
pixel size of 15 fim pixel -1 and a pixel scale of (X'13 pixel -1 . 
The 1.0" x 5.0" slit yielded a spectral resolution R ~ 24000 
(12 km s _1 ; the exact value depending on the binning used) at 
7000A, as measured from the full width at half maximum of 
the comparison lamp lines. Additional spectra were obtained 
over four consecutive nights in 2007 with the Inamori-Magellan 
Areal Camera & Spectr ograph (IMACS) M ulti-Object Echelle 
(MOE) spectrograph (Dress ier et al.l 1201 ll) on the 6.5 meter 
Magellan (Baade) telescope. The 0.6 " slit yielded a resolving 
power ofR = 21 000 (14 km s -1 ). 

The IMACS spectra were reduced and extracted with the 
Carnegie Observatories System for MultiObject Spectroscopy 
(COSMOS) pipeline (version 2.15), while the MIKE spectra 
with the MIK E python pipeline (version 2 .5.4), written by D. 
Kelson ( iKelson et alj 12000; iKelsor] [2003). The extracted or- 
ders for each star were averaged, normalized, and merged, 
yielding a wavelength range 5450A-9400A for the MIKE and 
6380A-8750Aforthe IMACS spectra, since t he large amount of 
reddening toward Westerlund 1 (Ay ~ 1 1 mag: lNegueruela et al.l 
l20T(ih limited the signal-to-noise ratio (S/N) at shorter wave- 
lengths. Heliocentric radial velocity corrections for each star 
were computed with the IRAfQ rvsao.bcvcorr routine and were 
taken into account in the subsequent radial velocity determina- 
tion. The spectra of Wdeb and W36 displayed double hydrogen 
and helium absorption lines consistent with a late-O or early- 
B spectral type, W13 displayed hydrogen and helium lines, in 
both absorption and emission, while WR77o displayed strong 
emission in He II XI6562, 8237 and N IV ^7103 -7128, consis- 
tent w ith the WN7o spectral type determined by ICrowther et al J 
(2006). Even though the S/N was not sufficient to detect the com- 
panion, we proceeded with the analysis of this single lined bi- 
nary. 



Table 1 : Log of Magellan spectroscopic observations 



UT Date Instru- Resolu- Disper- Total Exposure Time (sec) / S/N 



ment 


tion 


sion 










(A/AA) 


(A/pix) 




W36 


W13 WR77o 


01/03/2007 MIKE 


25500 


0.05 




1800/14 


1440/30 4800/7 


08/03/2007 MIKE 


24200 


0.10 


3600/22 


2700/38 


1440/50 4800/12 


09/03/2007 MIKE 


24200 


0.10 


3600/20 




1440/55 3300/14 


15/03/2007 MIKE 


24200 


0.10 




2700/42 


1260/49 3600/14 


28/06/2007 MOE 


21000 


0.25 






1800/31 ... 


29/06/2007 MOE 


21000 


0.25 




2400/20 


2400/35 . . . 


30/06/2007 MOE 


21000 


0.25 






3600/17 . . . 


01/07/2007 MOE 


21000 


0.25 




2400/35 


2400/7 


03/06/2008 MIKE 


24000 


0.10 


4800/13 


900/14 


1800/28 4800/7 


30/08/2008 MIKE 


23800 


0.15 




2700/56 


1440/81 3600/27 


01/09/2008 MIKE 


23800 


0.15 




2700/54 


1440/80 3600/20 


01/09/2011 MIKE 


24200 


0.10 


3600/32 







1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
NSF. 



Table Q] lists the UT dates of the observations, the spec- 
trograph and the corresponding resolution, dispersion and total 
exposure time and S/N for each target, determined at 7000 A. 
Typically three exposures of each target were obtained with the 
MIKE spectrograph to enable cosmic ray removal. The empty 
records indicate spectra with very low S/N, which made the de- 
tection of the lines and thus the use of those spectra impossible. 

3. Radial velocity analysis 

The radial velocities (RVs) were determined with ay 2 minimiza- 
tion method, which compares the observed spectrum to a grid of 
synthetic spectra. Each synthetic spectrum was constructed as 



Table 2: Radial velocity measurements 



Binary 


HJD 


Phase RVj ± a RVl 


(O-Ch 


RV 2 ± o- RV2 


(O - C) 2 




-2450000 




(km s~') 


(km s -1 ) 


(km s -1 ) 


(km s ') 




4167.87328 


0.42 


-188+35 


-19 


105+11 


- 16 




4168.86327 


0.64 


95±35 


56 


-165+11 


-17 




4620.60505 


0.22 


-205±35 


-40 


63±11 


-54 




5805.52158 


0.66 


130±35 


75 


-165+11 


-6 


W36 


4160.79354 


0.90 


100±32 


27 


-200+16 


-16 




4167.79030 


0.10 


-170±32 


-25 


110+16 


4 




4168.78891 


0.41 


-160±32 


-33 


95±16 


-10 




4174.78780 


0.30 


-200±32 


-3 


170+16 


-24 




4280.76847 


0.61 


100±32 


25 


-200+16 


1 




4282.66772 


0.21 


-180±32 


19 


165±16 


-34 




4708.56044 


0.09 






100+16 







4710.55897 


0.72 


100±32 


-27 


-270+16 


5 


W13 


4160.77158 


0.88 






-155+13 


7 




4167.76394 


0.64 


85±23 


-16 


-190+13 


-20 




4168.76923 


0.74 


145±23 


-2 


-210±13 


-1 




4174.76249 


0.39 


-180±23 


-7 


13+13 


-7 




4279.61631 


0.71 






-230±13 


-28 




4280.58040 


0.81 






-180+13 


17 




4281.62310 


0.92 






-125+13 


6 




4708.48711 


0.99 






-95+13 


-19 




4710.48815 


0.20 


-250±23 


-11 


75+13 


8 


WR77o 4160.84107 


0.98 






-25 ±30 


58 




4167.83559 


0.97 






-15 ±30 


97 




4168.82036 


0.25 






290 ±30 


-35 




4174.82619 


0.95 






-35 ±30 


108 




4282.49396 


0.54 






-230 ±30 


-100 




4620.66407 


0.61 






-340 ±30 


-74 




4708.52114 


0.57 






-250 ±30 


-61 




4710.51970 


0.14 






255 ±30 


13 



a sum of two models, each shifted over a range of radial ve- 
locities with a velocity step of 5 km s -1 . We used TLU STY 
models from the OSTAR2002 by lLanz & Hubenvl (l2003h and 
BSTAR2006 by lLanz & Hubenvl d2007l) grids for the compo- 
nents of each OB-type system. The parameters of the synthetic 
spectra for each system were selected to reproduce the observed 
spectra and the rotational velocities (v sin/) were adopted to be 
80 km s -1 and 1 10 km s -1 for the cooler and the hotter star of 
Wofs, respectively, and 110 km s -1 for the components of W36 
and W13. We inverted the TLUSTY model in order to match the 
emission line component of W13. Whilst this approach might 
seem somewhat unusual, it turned out to provide more stable 
results than a two Gaussian fit. The x 1 was computed as the 
quadratic sum of the differences between the observed and the 
synthetic spectrum for a range of radial velocities. We com- 
puted the minimization using 60A centered on each of the two 
available narrow helium lines (/LI6678, 7065), which yield more 
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accurate radial velocities compared to the broader hydrogen 
lines. For WR77o, w e used a WNE (early type) PoWR model 
dGrafener et al.l 120021; iHamann & Grafenerl |2003[ |2004|) . Each 
PoWR model depends on two parameters: the s tellar tempera- 
ture T * and the transformed radius R, (defined by Schmut z et al.l 
[1991 which takes into account the stellar radius, the mass-loss 
rate and the terminal-wind velocity. We adopted a model with 
a stellar temperature of T* = 50 100 K and transformed radius 
log R, = 0.7 and computed the minimization using 140 A cen- 
tered on the NIV AA7103 - 7128 line. The resulting velocities 
and errors for all four systems are given in Table |2] The fit pro - 
ducing the lowest x 1 yielded the best fit model of the grid, and 
therefore the radial velocities for each component. The errors 
were determined using the cr values, which were the mean values 
of the difference between the data and the model radial velocity 
curves determined in the following section. Figure Q] shows rep- 
resentative examples of the model fit to the observed spectra for 
each binary system, illustrating the good quality of the fit. 

4. Binary modeling 

With the radial velocity measurements at hand, we p roceeded 
to fit t he V, R and /-band light curves published by iBonanosI 
( 2007 ^ and the radial veloc ity curves using the PHOEBE inter- 
face (Pria & Zwitteil 12005 ) to the Wilson-Devinney code (WD; 
I Wilso n & Devinnev 11971 ). For all the binaries, we used the 
square root limb-darkening law, which works better for hot 
stars at optical wavelengths. Bolometric and passband limb- 
darkening coefficients were taken from the in-house computed 
limb-darkening PHOEBE 2011 tables. We fixed albedos and 
gravity brightening exponents to unity from theoretical values 
for stars with radiative envelopes. 

The errors for each p arameter were d etermined with the 
phoebe-scripter, following iBonan os (2009). The script uses the 
WD differential corrections algorithm to perform the minimiza- 
tion for the input values, adopts the minimizer results by apply- 
ing the corrections to the parameter table and repeats these steps 
1000 times. By calculating the x 2 minimization for each iter- 
ation we end up with statistics from which the errors for each 
parameter can be determined. 

Below we discuss each system separately and present the re- 
sulting orbital and astrophysical parameters of the components. 

4.1. W DEB 

The Westerlund 1 detached eclipsing binary (W^EfE!) is a 
double-lined, eccentric, detached system with a period of 4.4473 
d. Wdeb is located 4' south of the cluster core, but its colors are 
similar to those of the other cluster stars, a fact that indie ates that 
this binary is a member of Wd 1 . The components of Wdeb were 
the first main s equence stars to be identified in Westerlund 1 
dBona nos 2007j) and as such, their parameters can be compared 
to models of single star evolution and be used to estimate the age 
of the cluster. 

We obtained the first spectra for this system, which reveal 
an OB spectral type and confirm cluster membership. However, 
given the wavelength range available, a more accurate deter- 
mination of their exact spectral types and effective tempera- 
tures was not possible. We therefore adopted a temperature of 
25 000K in the analysis below as a starting point for our fitting. 
The few useful spectra (four) that we could work with, provide 
us with parameters that suffer from large errors (> 10%). 



Table 3: Results for Wdeb from combined LC and RV curve 
analysis 



Parameter 


Value 


Period, P 


4.4473 ± 0.0004 d 


Time of primary eclipse, HJD 


2453934.741 ± 0.003 


Inclination, i 


86 ± 2 deg 


Eccentricity, e 


0.175 + 0.003 


Longitude of periastron, u> 


252 ± 1 deg 


Surface potential, Cli 


5.9 ±0.1 


Surface potential, CI2 


6.50 ± 0.06 


Light ratio in V, L 2 /Li 


0.499 ± 0.005 


Light ratio in R, L2/L1 


0.501 ± 0.005 


Light ratio in /, L 2 /L l 


0.505 ± 0.005 


Mass ratio, q 


0.80 ± 0.07 


Systemic velocity, 7 


-40 ± 3 km s~> 


Semi-major axis, a 


34 ± 1 R 


Semi-amplitude, K l 


229 ± 19 km s" 1 


Semi-amplitude, K 2 


187 ±6 km s-' 


Radius 1 ", r l po i c 


0.202 ± 0.006 




0.208 ± 0.007 


r l.side 


0.203 ± 0.006 


fl.back 


0.206 ± 0.007 


-ri* 


0.204 ± 0.006 


Radius + , r^poie 


0.154 ±0.005 


r 2,point 


0.157 ± 0.006 


r2,side 


0.155 ±0.006 


r 2,back 


0.157 ±0.006 


r 2 * 


0.156 ±0.006 



Notes: 

f The radii are given in units of a. 
* Mean radius. 



We ran PHOEBE, using the three light curves and the ra- 
dial velocity curves, in the detached mode for the following free 
parameters: the inclination i, eccentricity e, the argument of pe- 
riastron to, the semi-major axis a, the radial velocity of the center 
of mass y, the surface potential for each star Q, the mass ratio 
<7=M2/Mi, the effective temperature T e ff, the period P, the time 
of primary eclipse HJDo and the band pass luminosity L2/Z4 for 
each band. We first fit the light curve to estimate P, HJDo, e, 
too, i, ^ for each star, then fixed these values and fit the radial 
velocity curve separately to estimate q, a and y. 

For a detached, eccentric system with partial eclipses such as 
Wdeb, the ratio of radii cannot be constrained only by photom- 
etry, so we used the spectroscopic light ratio to det ermine the 
primar y to be the larger and hotter star. Following North et al] 
d2010h . we next ran PHOEBE keeping the surface potential 
for the secondary component constant and then similarly deter- 
mined the surface potential for the primary component, while 
we adopted the value for q that was determined from the radial 
velocity curve. After the initial fit to the light curves and the ra- 
dial velocity curve, we iteratively refit the light curves and sub- 
sequently the radial velocity curve, until the values converged. 
The last step was to run PHOEBE for all the free parameters de- 
scribed above, but keeping q and Qi fixed. Table [3] presents the 
parameters resulting from the analysis with PHOEBE. The light 
and radial velocity curve fits are shown in Figure |2]and [3] respec- 
tively. The primary star, corresponding to the hotter component, 
was found to have a mass of 15 M Q and radius of 7 R Q while the 
secondary was found to have 12 M Q and 



RA=16:46:58.78, Dec=-45:54:31.9, J2000 
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I 1.0 E 

E 




0.7 ; 



WDEB 

. ' I... 
7050 



7060 7070 
Wavelength (A) 



7080 



1.0 



0.9 



0.7 




W36 

. i i i 
7050 



7060 7070 
Wavelength (A) 



7080 



1.6 - 




7045 7050 7055 7060 7065 7070 7075 7080 7060 7080 7100 7120 7140 

Wavelength (A) Wavelength (A) 



Fig. 1: Representative examples of the model fit (solid line) to the observed spectra (dotted lines), showing Hel 7065 for W13, 
W DEB , W36 and NIV 7103-7128 for WR77o. 




Fig. 2: Phased light curves of Wqeb in the I, R and V-band in units of normalized light, with best fit PHOEBE model overplotted 
(detached configuration). The apparent magnitudes outside of eclipse are / = 14.8, R = 17.35 and V = 20.3. The lower panel in 
each light curve shows the O - C residuals. 
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5 R Q . The \ogg values of the two stars were found to be ~ 4 
which confirms the main sequence nature of the components. In 
addition, the systemic velocity of y — -40 km s _I , which is simi- 
lar to that measured for the other three systems, provides further 
evidence for the membership of W deb to the cluster. Table [7] 
presents the physical parameters (mass, radius, logg, T e ff and 
log(L/L )) derived for each eclipsing binary. 




1 1 1 1 1 


1 1 1 1 1 1 1 


/i i\ 
nJ 1 \ 

\ / 


\ / \ I 


A 

'' l\I t/t 


1 1 A 















Fig. 4: Radial velocity curve for W36. The measurements are 
shown as filled circles for the primary and open circles for the 
secondary; overplotted is the best-fit model from PHOEBE, de- 
noted by a solid line for the primary and a dashed line for the 
secondary. 



Table 4: Results for W36 from combined LC and RV curve anal- 
ysis 



Fig. 3: Radial velocity curve for Wdeb- The measurements are 
shown as filled circles for the primary and open circles for the 
secondary; overplotted is the best-fit model from PHOEBE, de- 
noted by a solid line for the primary and a dashed line for the 
secondary. 



4.2. W36 

is the second brightest of our targets, with no previous 
spectroscopy available. It is a contact double-lined eclipsing bi- 
nary in a circular orbit with a period of 3.1811 d. The spectra 
of this system indicate that the components are of spectral type 
OB. As was the case for Wdeb, the determination of their exact 
spectral types and thus of their accurate temperatures was not 
possible. 

We ran PHOEBE, using the three light curves and the ra- 
dial velocity curves, in the contact mode for the following free 
parameters: the inclination i, the semi-major axis a, the radial 
velocity of the center of mass y, the surface potential for the pri- 
mary star £2i, the mass ratio q, the effective temperature of each 
star r e ff, the period P, the time of primary eclipse HJDo and 
the band pass luminosity Lq]L\ for each band. Table [4] presents 
the parameters resulting from the analysis with PHOEBE. The 
radial velocity and light curves are shown in Figure Hand re- 
spectively. 

The primary star of this system was found to have a mass of 
16 M Q and radius of 1 1 R G , while the secondary has 1 1 M Q and 
9 R Q (see Table [7). Since it is a contact system the masses that 
were determined are not the initial masses, as mass transfer has 
occurred. Values of log g < 4 found for both stars also testify to 
the evolved nature of the two stars. 



Parameter 


Value 


Period, P 


3.18110 ± 0.00003 d 


Time of primary eclipse, HJD 


2453909.8080 ± 0.0005 


Inclination, i 


73.0 ± 1.7 deg 


Surface potential, Sli 


3.29 ± 0.03 


Surface potential, Q.2 


3.28 (fixed) 


Light ratio in /, LalL\ 


0.514 ±0.005 


Light ratio in R, L2/Z4 


0.507 ± 0.005 


Light ratio in V, L 2 /L[ 


0.504 ± 0.005 


Mass ratio, q 


0.69 ± 0.03 


Systemic velocity, y 


-37 ± 2 km s- [ 


Semi-major axis, a 


27.5 ± 1.5 R a 


Semi-amplitude, K\ 


175 + 16 km s" 1 


Semi-amplitude, K 2 


235 ± 8 km s" 1 


Radius 7 , ri iPO i e 


0.378 ± 0.003 




0.47 ± 0.01 


reside 


0.398 ± 0.003 


r liback 


0.423 ± 0.004 


-ri* 


0.400 ± 0.004 


Radius 1 , r^poie 


0.316 ±0.003 




0.40 ± 0.01 


r2,sidc 


0.330 ± 0.004 


r 2,back 


0.357 ± 0.005 


r 2 * 


0.334 ± 0.005 



Notes: 

t The radii are given in units of a. 
* Mean radius. 



4.3. W13 

W130 is our brightest target: a near-contact double-lined eclips- 
ing binary in a circular orbit with a period of 9.266 d. The spectra 
indicate an emission plus absorption line object, consistent with 
the spectral types B0.5Ia+/WNVL and O9.5-B0.5I , respectively, 
report ed by the independent study of the system bv lRitchie et alj 
d2010h . The known spectral types provide additional constraints 



RA= 16:47:05.08, Dec=-45:50:55.1, J2000 



RA=16:47:06.46, Dec=-45:50:26.0, J2000 
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0.5 1 

I J I") ei s e 




0.5 1 

Ph ase 




0.5 1 

Phase 



Fig. 5: Phased light curves of W36 in the I,R and V-band in 
units of normalized light, with best fit PHOEBE model overplot- 
ted (contact configuration). The apparent magnitudes outside of 
eclipse are / = 13.1, R = 15.8 and V = 18.85. The lower panel 
in each light curve shows the O - C residuals. 



to determine the fundamental parameters of each component. 
In addition, the well determined temperatures of the members 
of W13 make this system the best available eclipsing binary in 
Wd 1 for a direct distance determination. 

We ran PHOEBE, using only the less noisy 7?-band light 
curve and the radial velocity curves, in the detached mode, as 
the components are not filling the Roche lobes, for the following 



free parameters: i, a, y, £2, q, P, HJDq, T e g 2 and L\ . Our radial ve- 
locity cu rves contain b o th our own RV measurements plus those 
made by iRitchie et al.l {2010). Table [5] presents the parameters 
resulting from the analysis with PHOEBE. 

Table 5: Results for W13 from combined LC and RV curve anal- 
ysis 



Parameter 


Value 


Period, P 


9.2665 ± 0.0003 d 


Time of primary eclipse, HJD 


2453902.3998 ± 0.0002 


Inclination, i 


61+2 deg 


Surface potential, Q.i 


4.6 ± 0.2 


Surface potential, Q.2 


5.6 + 0.2 


Light ratio in R, L 2 /L[ 


0.822 ± 0.005 


Mass ratio, q 


1.42 ±0.04 


Systemic velocity, y emis 


-48.2 km s- 1 (fixed) 


Systemic velocity, y a i, s 


-65.9 km s" 1 (fixed) 


Semi-major axis, a 


71 ± 1 R s 


Semi-amplitude, Ki 


198 + 7 km s-' 


Semi-amplitude, K% 


142 ± 4 km s" 1 


Radius 7 , ri iP0 i e 


0.31 ±0.01 




0.37 ± 0.04 


reside 


0.32 ± 0.02 


r l,back 


0.34 ± 0.02 


.ri* 


0.32 ± 0.02 


Radius 1 , r 2iPO ie 


0.29 ± 0.03 




0.31 ±0.03 


r 2.side 


0.30 ± 0.03 


r2,back 


0.31 ±0.03 


r 2 * 


0.30 ± 0.03 



Notes: 

t The radii are given in units of a. 
* Mean radius. 



The valu e s for the effective temperatures were taken from 
IRitchie et al.l d2010|) and correspond to the spectral types they 
determined. More specifically, T e ff l was fixed to 25 000 K and 
we let r e ff 2 be a free parameter with a starting value of 25 000 K. 
The light and radial velocity curves are shown in Figure |6]and [7] 
respectively. 

Moreover, for our analysis we adopted the y velocities for 
ea ch component ( y em = -48 km s _1 , y a b s = -66 km s ) as found 
bv IRitchie et al.l ((2010), and we transferred the whole system 
to velocity by adjusting the measured velocities according to 
the values above, since PHOEBE does not provide this option. 
In massive binaries, stellar winds can caus e the measured y to 
differ for each star (Ma ssev & Contilll977l) . This fact is highly 
correlated to outward atmospheric mass motion and given that 
the components are supergiants, it was important to include this 
effect. 

The component that is responsible for the emission in our 
spectra was determined to have 23 + 1 M , 23 + 2 R and 
a filling factor F = 0.96, values that are in agreement, al- 
th ough more precise th an the values independently determined 
bv IRitchie etal] J2010h : 23.2 +|-gM s , 22 + 2 R and F = 0.93. 
The component that causes the absorption in our spectra was de- 
termined to have 33 + 2 M H , 21 ±2Rp and F = 0.68, compared 
to the values found bv IRitchie etal] dioiOt) : 35.4 ^M , 21 + 
2 and F = 0.74. The high values of filling factors of the 
two components, make W13 a near-contact system. It is worth 
mentioning that our RV measurements alone yielded results very 
close to those above, with the star that causes the emission hav- 
ing a mass of 23 + 2 M G and radius of 22 + 2 R Q and the one 
that causes the absorption having 32 + 2 M and 23 + 2 R Q . 
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4.4. WR770 

WR77cQ is a single-lined eclipsing binary in a circular or- 
bit with a period of 3.52 d. The star that is visible in our 
spectra, is a Wolf-Rayet star, repo rted to have a spectral type 
WN6-7 dNegueruela & Clarkll2005l) . and later refined to WN7o 
( Crowthe ret al.ll2006l) . 

For a single-lined spectroscopic binary, one can only deter- 
mine the mass function of the system: 



f(M u M 2 ) = 



K\P 



2nG sin 3 i {M { + M 2 ) 2 



(1) 



where K 2 is the velocity semi-amplitude of the visible compo- 
nent. We first modelled the light curves with PHOEBE, for an 
estimation of the parameters i and P. Having determined these 
values, we estimated the mass function from K, P, G, i and found: 



Fig. 6: Phased /?-band light curve of W13 in units of normal- 
ized light, with best fit PHOEBE model overplotted (detached 
configuration). The apparent magnitude outside of eclipse is 
R = 14.65. The lower panel in the light curve shows the O - C 
residuals. 




Fig. 7: Radial velocity curve for W13. The measurements are 
shown as filled circles for the primary and open circles for the 
secondary; overplotted is the best-fit model from PHOEBE, de- 
noted by a solid line for the prima ry and a dash e d line for the 
secondary. The measurements from Ri tchie et al.l (1201 Oh for the 
primary and secondary are denoted by filled and open squares, 
respectively. All radial velocities have been shifted to a systemic 
velocity of zero. 



Furthermore we determined an inclination i — 61 + 2° , which is 
in good agreement with the v alue 62 ± 4° that w as determined 
from the binned light curve bv lRitchie et al.l(l2010h . The fact that 
we have almost doubled the number of RV measurements is the 
main reason that we determined a more accurate, but lower mass 
rat io q - 1.42 + 0.04 comparing to that of 1.53 + 0.1 determined 
bv lRitchie et alJ (l2010h . Finally, we improved the accuracy in the 
determination of the period to P — 9.266 5 + 0.0003 d compared 
to P = 9.2709 + 0.0015 d determined bv lRitchie etail (l2010h . 



(Mi + M 2 y 



= 21.0 ± 3.4 M Q 



(2) 



This large value implies a very large mass for the invisible 
component. If we further assume a mass of the Wolf R ayet star 
M 2 = 15.3 M Q as determined bv lCrowther et al.l (120061) for this 
system from the mass-luminosity relation of lSchaerer & Maederl 
(1 19921) for hydrogen-free Wolf-Rayet stars, we end up with a 
third degree equation, the solution of which gives us the result 
of a 40 M for the companion star. 

Table 6: Results for WR77o from combined LC and RV curve 
analysis 



Parameter 


Value 


Period, P 


3.520 ± 0.003 d 


Time of primary eclipse, HJD 


2453935.634 ± 0.003 


Inclination, i 


65 ± 3 deg 


Surface potential, Cli 


3.47 ± 0.03 


Surface potential, CI2 


2.78 + 0.03 


Light ratio in /, L2/Z4 


0.258 ± 0.005 


Light ratio in R, La/L\ 


0.250 ± 0.005 


Semi-amplitude, K 2 


360 ± 30 km s" 1 


Mass ratio, q 


0.37 (fixed) 


Systemic velocity, y 


-35 ± 3 km s~' 


Semi-major axis, a 


38 ± 2 R 


Radius 1 ", r lj)ole 


0.32 ± 0.02 


l"l,point 


0.33 ± 0.03 


reside 


0.32 + 0.02 


Ij.back 


0.33 ± 0.03 


-ri* 


0.32 ± 0.03 


Radius f , r 2>po i c 


0.24 ± 0.02 


r2,point 


0.28 + 0.03 


r 2,sidc 


0.25 ± 0.01 


r2,back 


0.27 ± 0.02 


r 2 * 


0.25 + 0.02 



Notes: 

f The radii are given in units of a. 
* Mean radius. 



We can continue with this preliminary analysis to constrain 
some of the parameters and find a more rob ust solution using 
PHOEBE. Using the following values from Crowther (2007): 
log(L/L ) = 5.54 and T eS = 50000 K for a WN7 spectral type 
Wolf-Rayet star, we find the radius of the star to be ~ 8 R Q . 



RA=16:47:05.37, Dec=-45:5 1:04.7, 12000 
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We ran PHOEBE using the / and /?-band light curves and 
the radial velocity curve of the secondary component, in the de- 
tached mode, for the following free parameters: i, a, y, Q., T e g, 
P, HJDq and L\ . The eccentricity was fixed to and the mass 
ratio q to 0.37, as determined above. 




0.5 1 

Phase 



Table [6] presents the parameters resulting from the analysis 
with PHOEBE, which yielded a stable solution in agreement 
with our original estimates. The light and radial velocity curves 
are shown in Figure [8] 

The Wolf-Rayet star was found to have 16M and 9.7 R Q 
while the primary companion 43 M and 12.3 R e . The T e ff of the 
Wolf-Rayet star was found to be 44 000 K which is quite close 
to the expected one for its type. Even though the companion is 
characterized by a higher T e ff (64 000 K) and a larger mass and 
radius, it is not visible in our spectra due to the low S/N, which 
adversely affects absorption lines. The high values of the filling 
factors, F = 0.75 and F = 0.85, respectively, indicate a post- 
mass transfer system. 

Wolf-Rayet stars are intrinsically variable, therefore part 
of the observed scatter in the light curve is due to this effect. 
Mod eling, which also t akes into account atmospheric effects 
(e.g. iPerrier et al.l |2009), would possibly improve the fit and 
yield more accurate parameters. Furthermore, a more precis e de- 
termination of the period would likely correct the slight phase 
shift of the radial velocity curve with respect to the model. We 
conclude that due to all the assumptions made for WR77o, the 
estimated parameters should be used with caution. 

Table 7: Physical parameters of the four Wd 1 eclipsing binaries 




Fig. 8: Phased light curves of WR77o in the / and /?-band 
in units of normalized light and radial velocity curve, with 
PHOEBE models overplotted (detached mode). The apparent 
magnitudes outside of eclipse are / = 14.3 and R - 17.55. The 
lower panel in each light curve shows the O - C residuals. 



Parameter 




W36 


W13 


WR77o 


Mi (M ) 


14.8 ±3.5 


16.3 ± 1.5 


23.1 ± 1.1 


43.4 


± 6.8 


M 2 (M Q ) 


11.9 ± 1.8 


11.3 ± 1.8 


32.9 ± 1.9 


16.1 


±2.5 


Ri (Ro) 


6.9 ± 2.0 


11.0 ± 1.2 


23.0 ± 1.5 


12.3 


± 2.0 


R 2 (Ro) 


5.3+2.0 


9.2 ± 1.2 


21.3 ± 1.5 


9.7 


t 2.0 


loggi 


3.93 ± 0.27 


3.57 ±0.10 


3.08 ± 0.06 


3.89 


t0.16 


l0gg 2 


4.07 ± 0.33 


3.56 ±0.13 


3.30 ± 0.05 


3.67 


t0.19 




29 000 ± 580 


30 000 ± 1 500 


25 000 (fixed) 


64000 


±5 000 


7-eff 2 (K) 


27 000 ± 570 


25 500 ±2 900 


25 000 ± 1 000 


44000 


±4500 


logLi/Lo 


4.48 ± 0.25 


5.19 ±0.13 


5.27 ± 0.06 


6.36 


t0.20 


logL 2 /L Q 


4.13+0.33 


4.51 ±0.23 


5.21 ± 0.09 


5.50 


t0.25 



5. Evolutionary status 

For a comparison of our measu red parameters with ev olutionary 
models, we used the models of Ekstrom et al ] (l2012h . for single 
stars in the range 0.8 - 120 M at solar metallicity (Z=0.014), 
which include rotation. These evolutionary models are not ap- 
propriate for post mass-transfer binaries, however we proceed 
with the comparison, as it provides useful information, while 
computing binary evolution models is outside the scope of this 
work. Figure [9] compares the parameters of each binary (masses 
vs. radii) with evolutionary tracks and isochrones for single ro- 
tating stars at solar metallicity, while Figure [10]plots them on an 
H-R diagram. 

The position of the components of W36 on the mass-radius 
diagram (Figure [9} indicates an age of > 10 Myr for the sec- 
ondary, while the discrepancy is smaller in the H-R diagram. The 
appa rent conflict with pre vious studies that find an age of 4.5-5 
Myr dCrowther et al.ll2006l) can be resolved by considering that 
W36 is a contact binary, which has undergone mass transfer. On 
the other hand, ~Wdeb with its unevolved members of 15 M and 
12 M is, in principle, ideal for estimating the age of the clus- 
ter. Its position on this diagram supports an age that is consistent 
within errors to previous studies, as detached unevolved systems 
can be treated as single stars. The position of the secondary of 
W 1 3 also indicates a similar age for Wd 1 . WR77o is not plotted 
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Fig. 9: Comparison of the parameters of W deb (stars), W36 (cir- 
cles) and W13 (squares) with evolutionary tracks (solid lines) 
for stars with 5-60 M and isochrones (dotted lines) for sin - 
gle rotating stars at solar metallicity from lEkstrom et afl (120121) . 
Isochrones, from the bottom up, correspond to 3, 4, 5, 6, 8 and 
10 Myr. Single star isochrones are not compatible with the mea- 
sured parameters for W36 and W13, indicating mass transfer. 



in the diagram since we were unable to detect both stars spectro- 
scopically to determine accurate parameters. 

The position of the components of W36 on the H-R diagram 
(Figure [TOl show that both stars appear overluminous for their 
mass, which provides further evidence th at the system has un- 
dergone mass transfer. [North et all (1201 Ol) have also reported in- 
stances of non-detached eclipsing systems whose components 
appear overluminous or underluminous for their mass when 
compared with tracks for single stars. The position of Wmj is 
in good agreement with the single star models although the cur- 
rent error bars do not allow a precise age estimation. Finally, the 
primary of Wl 3 is in agreement with the models, while the sec- 
ondary appears underluminous, once again providing evidence 
for mass transfer. 



Fig. 10: Comparison of the parameters of Wdeb (stars), W36 
(circles) and W 13 (squares) with evolutionary tracks (solid lines) 
and isochrones (dotted lines) for sing le rotating stars at solar 
metallicity from lEkstrom et au (12012b . Isochrones, from left to 
right, correspond to 3, 4, 5, 6, 8 and 10 Myr. The comparison 
finds Wqeb in good agreement with single star models, the com- 
ponents of W36 to appear overluminous for their mass, and the 
primary of W13 in agreement, while the secondary appears un- 
derluminous. 



direct, independent eclipsing binary distance measurements to 
Wd 1. 

Of the four systems, W13 has the most accurately deter- 
mined radii and effective temperatures, and so we proceed to 
comp ute the distance to W esterlund 1 based on this eclipsing bi- 
nary. [Ritchieletal] (|20 10) determined accurate spectral types for 
the components of W13: B0. 5 Ia + /WNVL for the emission line 
component and O9.5-B0.5I for the absorption line component, 
which correspond to a T e ff = 25 000 ± 2 000 K. We adopted this 
conservative error for both components for the distance estima- 
tion. Using the radii resulting from our analysis, we compute 
log(Li/L ) = 5.27 + 0.20 and log(L 2 /L ) = 5.20 + 0.20 , and for 
the combined light: log(L tota i/L ) = 5.54 ± 0.11. The absolute 
magnitude of a star at a certain wavelength A is given by: 



6. Distance to Westerlund 1 



The distance to Westerlund 1 has previously been estimated to 
be ~ 4 - 5 kpc, using a variety of methods, e.g.: WR stars 
(5.0 + °; g kpc:ICrowtheretaill2006l) . HI observations (3.9 + 0.7 
kpc; Kothes & Dougherty 2007), OB supergiants (> 5 kpc; 
Negueruel a et alJl2010h. a nd pre-main sequence stars (4.0 ± 0.2 
kpc; [Gennaro et al.ll20 1 lb . Eclipsing binaries provide powerful 
tools for measurin g independent an d accurate distances to stel- 
lar clusters (e.g. |Rauw et aD 120071) and nearby galaxies (e.g. 
Bon anos et al.l 120061 1201 llh The analysis of the four systems 
described in the previous sections therefore opens the way for 



M, =M, 



bol 



O 



-BC A - 2.5 log — 



(3) 



where M Q bo1 = 4.75 mag is the bolometric luminosity of the Sun 
and BCa t he bolometric correc tion. We used theoretical BC val- 
ues from lMartins et al.l (120061) . whi ch are based on state-o f-the- 
art CMFGEN model atmospheres dHillier & MiUerlll998l) . For 
an 09.51 star, the theoretical BCy = -2.62 mag, which yields 
My = -6.47 mag. 

The optical out-of-eclipse magnitudes of W 1 3 are V = 
17.50, R = 14.66 and / = 12.02 dBonanosll2007l) IClark et al.l 
d2010 ) found the V-band calibration of iBonanosI (120071) to be 
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offset to that of lClark et all (120051) . while the R and /-band cal- 
ibrations were in good agreement, and reported magnitudes for 
W13 of B - 21.1, V = 17.19,/? = 14.63,/ = 12.06, corre- 
sponding to phase ~ 0.25. For the distance estimation , we adopt 
the higher-quality photometry of IClark et ail (120051) . obtained 
with a larger telescope and a conservative error of 0.04 mag 
in each band, which corresponds to the difference in the R and 
/-band measurements. 

The distance d is estimated from the distance modulus equa- 
tion: 



My 



5 x log 



\l0pc 



(4) 



where m v is the apparent V-band magnitude, and Ay = Ry X 
E(B - V) is the extinction. From the photometr y of IClark et al.l 
d2005l) and adopting (B-V) = -0.26 from lMartins et al.l (120061) . 
we obtain E(B - V) = 4.17 + 0.05 mag, which yields Ay = 
12.93 +2.09 mag (forRy = 3.1 +0.5), in agre ement with the av- 
erage value found of A y = 13.6 mag found bv lClark et al] d2005l) 
for OB supergiants. Equation |4]therefore yields d = 1.4^g kpc, 
which is a very imprecise estimate due to the large uncertainty 
in Ay, demonstrating that near-infrared photometry and an ac- 
curate determination of the extinction is essential for an accurate 
distance measurement. 

We therefore proceed to use near-infrared photometry, 
which is less affected by dust extinction. The available 2MASS 
hotometry of W13 (at an unknown phase) is J — 9.051 +0.020 
and we adopt a conservative error of 0.16 mag (the depth of 
W13 primary eclipse), i.e. J = 9.05 + 0.16 mag. 

For an 09.51 star, the theoretical BCy= -3.24 mag, which 
yields Mj = -5.85 mag. We use the following equation to esti- 
mate the distance: 



mj - Mj - Aj = 5 x log 



/ d 
\ \0pc 



(5) 



We adopt Aj/A K = 2.50 + 0.15 (llndebetouw et al.ll2005l) and 
the average value for Ak = 0.82 mag, determined toward the 
W R star "R" (or WR77q) , a close neighbor of W13 ~ 6" away 
by ICrowther etail (120061) . We note that Star "U" (or WR77s), 
also nearby, has Ak = 0.74 mag, therefore we adopt a conser- 
vative error in Ak of 0.1 mag. We find Aj — 2.05 + 0.28 mag 
and derive a distance: d = 3.71 + 0.55 kpc to W13, which is 
the first eclipsing binary distance to the cluster. This value with 
an accuracy of 15% is in good agreement with recent estimates, 
which is noteworthy as eclipsing binaries are reliable, direct and 
independent distance indicators. For future determinations of the 
true distance modulus (DMq) with better estimates of the extinc- 
tion, we also report the value DMq + Aj = 14.90 + 0. 16mag. The 
precision of our measurement can be improved to < 5% by: (a) 
future near-infrared photometry to unambiguously disentangle 
the extinction and reddening law toward W13, and (b) spectral 
energy distribution (SED) modeling to improve the temperature 
determination. 



7. Discussion and conclusions 

We have presented accurate (better than 15%) fundamental pa- 
rameters of four massive eclipsing binaries in the young mas- 
sive cluster Westerlund 1, confirming they are all cluster mem- 
bers and demonstrating once more the success of the survey 



strategy. The mass range found for the eight component stars 
is 1 1 -43M , and three of the systems (W13, W36 and WR77o) 
seem to be post mass-transfer binaries. Our results for W13 are 
in agreement with t hose presented in the independent study by 
Ritchi e et al. (2010). We were able to improve the accuracy of 
the parameters of W13 and therefore provide confirmation on 
the dynamical constraint on the high proge nitor mass (> 40 M ) 
of the ma gnetar CSO J 1 647 10 .2-4552 16 dMuno et alj|2006l) re- 
ported by Ritc hie et al.l (HoiO), which accounted for binary evo- 
lution. Such massive stars are believed to form black holes rather 
than neutron stars, at least in a single star scenario. However, 
a strong interaction, as that of mass stripping in a binary sys- 
tem, can change the progenitor mass range significantly, without 
causing any problem in o ur understanding of stellar evolution 
dBelczvnski & Taam 2008). Our discovery of a massive ~ 40 M 
component for WR77o is also significant as it will serve as a 
second dynamical constraint on the mass of the progenitor of 
the magnetar. Future high S/N spectra should reveal the massive 
component, which, combined with a new, higher quality light 
curve, will enable an accurate determination of the component 
masses. 




10 

Mass (M Q ) 



Fig. 11: Mass and radius determinations of stars in eclipsing bi- 
naries, accurate to < 10% and co mplete > 3OM compiled from 
the literature by iBonaiiosI (120091) . now updated with the massive 
secondary of W13. The s olid line is the Z _ .02 zero-age main 
sequence (ZAMS) from Scha ller et all dl992h ; the dashed line is 



the Z = 0.008 ZAMS from Schae rer et al.ld!993b : the dotted line 



The H and K bands only have upper limits. 



is the terminal age main sequence (TAMS). W13 brings the total 
number of stars fulfilling the requirements of mass and accuracy 
to 18. 



W13 is also noteworthy because the a b sorption compo nent 
(33 + 2M Q ) fulfills the criteria set forth bv lBona nos (2009| and 
increases the total number of very massive stars (> 30 M Q ) with 
< 10% measurements from eclipsing binaries to 18. Figure [TT] 
shows the updated mass radius diagram with eclipsing binary 
measurements accurate to 10% and complete > 3OM . The slow 
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rate of increase of this sample (one star per year) is a reflection 
of the difficulty of identifying very massive stars and obtain- 
ing follow-up observations. The availability of eclipsing bina- 
ries with accurate parameters now enables direct, independent 
eclipsing binary distance determinations to the cluster. We have 
used our results for W 1 3 to estimate the first eclipsing binary dis- 
tance to the cluster and therefore the magnetar, finding a value 
of d — 3.7 + 0.6 kpc, in agreement with previous determinations 
using other methods. Future near-infrared photometry, a better 
light curve, in addition to careful modelling of the SED has the 
potential to provide a measurement accurate to < 5%, which will 
be crucial for our comprehension of magnetars. 

Lastly, we find Wdeb, despite its distance from the core of the 
cluster, to have the same systemic velocity as the other eclips- 
ing systems, indicating it was ejected due south. Future proper 
motion studies will determine the time of ejection and possible 
connection with other massive cluster stars located outside the 
core, such as the ma gnetar, luminous blue v ariable, WR77aa (T) 
and WR77b (N; see ICrowther et al.ll2006h . possibly indicating 
frequent ejections of stars from the cluster c ore, as was reported 
for Westerlund 2 dRoman-Lopes et aLll201 ll) . 

In conclusion, we note that young massive clusters remain 
largely unexplored regarding their binary star content, although 
binaries are extremely useful tools for providing constraints for 
models of evolution and formation for single and binary stars, 
at a range of metallicities. In particular for Wd 1, we have used 
the four known eclipsing binaries to obtain fundamental param- 
eters for eight massive stars and to derive the first independent 
distance determination to the cluster. 
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